Robustness of the helical edge states in topological insulators 
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Topological insulators (TI) are materials having an energy band gap in the bulk and conducting 
helical electronic states on the surface. The helical states are protected by time reversal symmetry 
thus are expected to be robust against static disorder scattering. In this work we report atomistic 
first principles analysis of disorder scattering in two-probe transport junctions made of three dimen- 
sional TI material Bi2Se3. The robustness of the device against disorder scattering is determined 
quantitatively. Examining many different scattering configurations, a general trend emerges on how 
strong the perturbing potential and how it is spatially distributed that can derail the helical states 
on the Bi2Se3 surfaces. 

PACS numbers: 73.43.-f, 72.25.Hg, 72.10.Fk, 73.20.-r 
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The concept of topological insulator (TI) has attracted 
tremendous recent attention (Tj-Q due to its fundamental 
role in classifying electronic structures of materials hav- 
ing a band gap. TI is also important for potential appli- 
cations in nanoelectronics due to its peculiar transport 
properties. TI has an energy gap in its bulk band struc- 
ture but has metallic helical states on its surface^, 0]- 
Importantly and interestingly, the direction of electron 
spin of the helical states is locked to perpendicular to the 
electron momentum k: electrons moving in the positive 
direction of k have their spins pointing to one direction 
while those moving in negative k pointing to exactly the 
opposite direction. The helical states in TI is induced by 
spin-orbit interaction (SOI) and are protected by time- 
reversal symmetryplQ, as a result electron back scatter- 
ing by disorder is suppressed because the electron's spin 
cannot be flipped if the disorder does not break time- 
reversal symmetry. The robustness of the helical con- 
ducting channels against disorder scattering is a most 
distinct characteristic of TI and is the most important 
factor for practical applications of TI to nanoelectronics. 

An important question therefore arises: how robust are 
the helical states in real TI materials? In other words, 
while weak disorder is not expected to back scatter elec- 
trons in the helical states, a strong disorder may mix the 
helical states with the bulk states of the disordered ma- 
terial and prevent helical states from establishing in the 
first place. Understanding how electrons traverse disor- 
dered region is urgently needed for establishing a phys- 
ical picture about TI based nanoelectronics. This is es- 
pecially important in light of the recent experiments in 
two-probe TI transport junctions in which microscopic 
transport parameters were measured[6J|. To the best of 
our knowledge, we are not aware of previous analysis on 
the robustness of helical states from atomic first princi- 
ples that fully takes the microscopic electronic structures 
of the real material into account. Here, by examining 
many different scattering configurations ab initio, a gen- 
eral trend is discovered on how strong a perturbing po- 



tential and how it is spatially distributed that can derail 
the helical states. 

To be specific, we consider the popular three- 
dimensional (3D) TI material Bi2Se3 forming a two-probe 
transport junction in the form of a two-dimensional (2D) 
film as shown in Fig[T] Transport is in the x-direction, 
the film is periodic in the y-direction, and the thickness 
of the film (z-direction) is six quintuples or 30 atomic 
layers (AL) which is thick enough to isolate interactions 
between the top and bottom surfaces across the fllm0- 
9]. In momentum space (k-space), the x-direction cor- 
responds to the momentum connecting the T-point to 
the M-point (FigJTJ:). The two-probe device consists of 
three parts: left/right electrodes and the scattering re- 
gion. The scattering region (indicated by letter "C" ) has 
two primitive cells (PC) along the x-direction, a total of 
60 atoms[10]. Potential perturbation SV mimicking dis- 
order is added in the scattering region of the device. Both 
left /right (L/R) electrodes are perfect Bi2Se3 films ex- 
tending to x = ±oo. If there is no SV, the two-probe de- 
vice reduces to a perfect infinitely big Bi2Se3 film. With 
SV, we investigate the robustness of the incoming helical 
state as it traverses through the scattering region. 

Our quantum transport calculation is based on first 
principles where density functional theory (DFT) is car- 
ried out within the nonequilibrium Green's function 
(NEGF) formalis m [111 , as implemented in the Nanodcal 
transport package |12l fl3l |. Very briefly, in the NEGF- 
DFT technique, the Hamiltonian of the open device 
structure is determined by DFT, the density matrix 
which includes nonequilibrium quantum statistical infor- 
mation of quantum transport is determined by NEGF, 
the transport and electro-static boundary conditions of 
the device are handled by real space numerical tech- 
niques. For further details we refer interested readers 
to the original liter atureflll]. 

After the two-probe device Hamiltonian of a perfect 
junction (Figfjji) is converged, we analyze disorder scat- 
tering using a diagonal disorder model where a perturb- 



2 



L C R i z 




C2 C3-C5 



FIG. 1: (color online) (a) Atomic structure (one period along 
y-axis) of the two-probe Bi2Se3 transport junction. Boxes 
labeled L, C, and R form the two-probe unit cell composed 
of the left electrode, scattering region, and right electrode. 
Atoms to the left of L and right of R form device electrodes 
which extend to x = ±00. The thickness of the Bi2Se3 2D 
film is 6 quintuples terminating with Se atom. Each of the L, 
C, R boxes has two primitive cells (PC). The framed area in 
green between the first and second atomic layer (from bottom) 
is the area where electron wave function is shown in Fig[3] 
The disorder configurations are given in Table [I] Region C for 
configurations C2, C3-C5 is shown in lower left panels where 
SV is added to the atoms in green, (b) Top view of the lattice 
vectors [lOJ ai and a 2 and the atom positions Ai, A2, and A3. 
The dotted rectangle shows the size of each region in the x-y 
plane, (c) The Brillouin zone in the k x -k y plane at k z = 
and the special k points P, M, and K. 



ing potential SV is added to the diagonal elements of the 
Hamiltonian matrix corresponding to the scattering re- 
gion (region C in FigOJi). This approach is adequate for 
our purpose while avoiding the prohibitively large com- 
putation if SV were to be treated self-consistently. In the 
study, we investigate many SV configurations as listed in 
Table-IH 

Fig UK plots the calculated band structure of the CI 
configuration (no disorder) - only the energy along the 
T-M direction (see FigOJ,) is shown. The dotted lines 
are bulk bands and the solid line is the Dirac band of 
the helical states of the surfaces. The calculated Fermi 
level (vertical line) locates at round 0.072eV above the 
Dirac point. A gap of ~ 0.27eV separates the bulk va- 
lence and conduction bands while the Fermi level cuts 
the Dirac band. These agree with previous ab initio 
results 0- For CI, the Fermi wave vector is calculated to 
be 0.0333/A( see Figl2k). The measurements on Bi 2 Te 3 
gave the same value [6|. 

For the 2D film, the transmission coefficient T(E, k y ) 
is a function of the incoming electron energy E and its 



Label 


SV(eV) 


Configuration 


CI 





perfect junction 


C2 


0.01 


15 atoms, lower half PC (1-15 ALs) 


C3 


0.01 


30 atoms in PC (1-30 ALs) 


C4 


0.1 


30 atoms in PC (1-30 ALs) 


C5 


1.0 


30 atoms m PC (1-30 ALs) 


C6 


[0,0.1] 


60 atoms in C (1-30 ALs) 


C7 


[0,0.5] 


60 atoms in C (1-30 ALs) 


C8 


[0,U] 


60 atoms in C (1-30 ALs) 


C9 


[0,1] 


atoms in any single quintuple 



TABLE I: Disorder configurations investigated. Configura- 
tion means where the diagonal perturbing potential SV is 
added. In C2 to C5, the same value of SV is used for all 
atoms specified in the configuration. PC is the first primitive 
cell in region C (C has two primitive cells, Fig[T^). In C6 
to C9, the value of SV applied independently for each atom 
is drawn randomly with equal probability from the energy 
ranges shown. 

transverse momentum k y , and a 2D transmission T(E) is 
obtained by integrating T{E 1 k y ) over k y . We shall first 
investigate T(E, 0) which is the transmission for electrons 
incoming along the T-M direction (k y — 0). We found 
that T(E, 0) has an even-integer step structure shown by 
the solid line in FigJ2t>,c. This is because there is no scat- 
tering in CI and every incoming Bloch wave (at energy 
E) perfectly transmits through, hence T(E, 0) must be 
an even integer (even due to spin degeneracy). In partic- 
ular, between the bulk gap, T(E, 0) = 2 for all energy E 
because there is only the doubly degenerate helical state 
in that energy range. 

Extraordinarily, T(E, 0) is still 2 in the bulk gap even 
if there is a moderate disorder in the scattering region: 
as shown for C3 (dashed), C4 (dotted), and C6 (dash- 
dotted) configurations in FigJ5}3. The results of C2 is very 
close to that of C3 thus not plotted. These indicate that 
SV cannot cause back scattering to the incoming helical 
states and transport is immune to it. In contrast, out- 
side the bulk gap, i.e. for E < -0.009eV or E > 0.27eV 
where bulk bands contribute to transport, SV reduces 
the total transmission. Most surprisingly, even for the 
strong scattering C5 (dotted in FigJU;), T(E,0) = 2 for 
the range E = [0.1, 0.27]eV. Here, SV does reduce trans- 
mission at lower energies E < O.leV, indicating that a 
strong perturbing potential can prevent the surface he- 
lical state from establishing. Our calculation also shows 
that the helical states are stable against a SV as strong 
as leV if it is added to all atoms in any single quintu- 
ple of a PC in region C (configuration C9). However, a 
SV = 0.3eV added to all atoms in two quintuples of a 
PC in region C might collapse the helical states. This 
indicates that having disorder to a larger volume is more 
detrimental to the helical states than increasing the dis- 
order strength in a smaller volume. For strong random 
potentials of C7, the transmission (dashed line in FigJSJ:) 
is only slightly reduced from 2 at lower energy but the 
helical state is destroyed near the Fermi energy. These 
clearly demonstrate the robustness of the helical states 
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FIG. 3: (color online) The real-space scattering states WF1 
(left panels) and WF2 (right panels) in regions L, C, and R of 
the device for configurations CI, C5, C6, and C7 (from upper 
to lower panels). 



FIG. 2: (color online) (a) Energy bands of the helical states 
(solid) and the bulk states (dotted) in the Bi2Se3 film between 
the k points M and F. The Fermi level is indicated by the ver- 
tical line, (b) Transmission T(E, 0) versus electron energy E 
near the Fermi energy Ef for electrons transporting along the 
F-M direction in configuration CI (thin solid), C3 (dashed), 
C4 (dotted), and C6 (dash-dotted). Inset: 2D conductance G 
versus the electron energy in unit of the conductance quanta 
Go = %- in configurations CI, C4, C6, C7, and C5. (lines 
from up to down), (c) T(E,0) versus E for strong disorder 
configurations C5 (dotted) and C7 (dashed). Inset: transmis- 
sion at Fermi energy E=0 (dotted, blue) and E=0.1 eV (solid, 
red), versus the strength U of random potential (C8). 

against disorder scattering. 

An important character of helical states for perfect TI 
(CI) is its spin-momentum locking, namely the deviation 
angle a of the electron spin u from perpendicular to its 
direction of momentum vanishes to a high precision [7] : 
a = er ■ = 0. In this work, we compute er by spatially 
averaging spin in the L-C-R region (see FigO} using the 
method described in Ref.[0j], and fix k along the trans- 
port direction x. With disorder, we found the angle a 
may become nonzero. In C5 and C7 where the helical 
states are destroyed by strong disorder, a is found to 
be as large as 0.1, ten times greater than that of other 
configurations. 

To gain further insight, let's consider wave functions of 
the helical states after they suffer the SV scattering. In 
perfect Bi2Se3 (CI), there are various Bloch states and 
also localized evanescent states at any given energy E. In 
the bulk gap, the Bloch states are two degenerate helical 
states that live on the bottom (wave function denoted as 
WF1) and top (WF2) surfaces, respectively, and WF2 
has opposite average spin compared to WF1. During 
transport, the incoming helical states must traverse the 
scattering region (see FigQJi). In disorder configurations 
(C2-C9), incoming electrons in a Bloch state of spin a 



in the nth channel 4>f = 4>f a (subscript I indicates left 
incoming), are scattered by SV. The scattering process 
couple the incoming wave to the evanescent states and 
other propagating channels, the electrons are reflected or 
transmitted to the outgoing modes of index m, i.e. ipf 1 
or tp™, respectively. The scattering wave function ^ 
corresponding to the n-th incoming channel is calculated 
as[nj]: 

f # + Em CP?* inL 
n = { r c in C (1) 

where ^ is expanded in terms of the ex- 
tended/evanescent modes m. After the Hamiltonian of 
the two-probe device is converged by the NEGF-DFT 
self-consistent calculation, it is acted on ^ for a given 
energy E, and the resulting algebraic equations for the 
transmission (ij^) and reflection (r^J amplitudes are 
calculated so that ^ can be determined [111]. 

FigO plots the real space WF1 (left panels) and WF2 
(right panels) sliced in the x-y plane between the first 
and second ALs near the bottom surface (green frame in 
FigHt-)- The electrons are concentrated around Ai (Se 
atom on the first AL) and A2 (Bi atom on the second 
AL) positions. The first row shows results for the per- 
fect film CI. On average WF1 has an amplitude two 
orders of magnitude larger than WF2 since WF1 (WF2) 
is confined to the bottom (top) surface. In C2 which has 
relatively weak disorder on the 15 atoms in the lower half 
of the PC, WF1 is affected but WF2 remains the same 
as that of the CI - because WF2 lives on the top sur- 
face away from the atoms perturbed by SV . For C3 and 
C4 where SV is added across the film, both WF1 and 
WF2 are significantly affected in region C but much less 
so out side it, hence the transmission and the 2D con- 
ductance is not affected by SV (Figj2j) . The patterns of 
real space WF1 and WF2 for C2-C4 are similar to those 
of CI (not shown). For C5 (second row in Fig|3|) which 
is a strong SV case, WF1 and WF2 are no longer sepa- 
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rated in real space and have similar orbital distributions 
with opposite spin orientations. The two states are now 
mixed through the him with almost symmetric profiles 
of amplitude along z between the two surfaces. Namely, 
scattering states in regions L, C, R are all significantly 
affected thus transmission reduces (indicated in Fig 12]). 
Our calculation indicates that a SV of ~ 0.3eV destroys 
the helical states if being applied to a region enclosing 
more than ~ 10 atoms. If SV is added to any single 
atom - a single impurity, the helical states arc robust up 
to SV ks leV which is a very large local perturbation in- 
deed. Since there are 60 atoms in region C, 10 perturbed 
atoms correspond to an impurity concentration of 16.7%, 
while a single impurity corresponds to a concentration of 
1.67%. 

Similar robustness of the helical states is found for ran- 
dom potentials: for C6 (third row), both WF1 and WF2 
are modified from the perfect film case but they are still 
localized in their respective surfaces thus transmission 
T = 2 in the bulk gap is largely intact. For C6 the aver- 
age perturbing potential SV = 0.05eV, not strong enough 
to derail the helical states. For C7 - stronger random po- 
tential (fourth row), WF1 and WF2 become similar in 
real space (i.e. not confined in surfaces) indicating the 
helical states are destroyed. It is very interesting to note 
that even though SV is added only in region C, the scat- 
tering states (WF1, WF2) are affected in all regions L, 
C, and R reflecting the scattering processes. 

To correlate the robustness of transmission T(E, 0) 
against the strength of the disorder potential SV, we 
have calculated configuration C8 (see TablelT]) where 
SV is randomly drawn from ranges [0, U], results are 
shown in the inset of Fig(2£ versus U. At Fermi energy 
E = (blue dotted line), T(E,Q) = 2 until U reaches 
~ 0.2eV at which T(E, 0) substantially reduces (inset 
of Fig[2}:). By investigating which atomic orbitals con- 
tribute to the scattering wave functions and also the real 
space WF1/WF2, the scattering states are found to be 
unconfined from the surfaces at this value of U, namely 
the helical states are destroyed. When U increases fur- 
ther, T(E, 0) oscillates by the random disorder scatter- 
ing configurations. At higher energies, we found that 
the helical edge states are much more robust. For in- 
stance, at E — O.leV (solid curve in the inset of FigJ5b), 



T(E,0) ~ 2 until U ~ 0.5eV. This is reasonable since 
electrons with higher energy are less affected by poten- 
tial perturbations. 

So far we focused on transmission in the x-direction, 
T(E,k y = 0). For the film, a 2D transmission T(E) is 
obtained by integrating T(E, k y ) over k y , which gives the 
2D conductance G(E) = T{E) x G a where G a = e 2 /h. 
The inset of FigJ^Jj plots G(E). For energies outside 
the bulk gap, e.g. E < or E > 0.27eV, bulk propa- 
gating modes participate transport and therefore G de- 
creases according to the increase of the disorder. Within 
the gap, G{E) changes negligibly for weak disorder C4 
and C6, but more significantly in strong disorder C5 and 
C7. In contrast to T(E,0) (the main figure in Figl5]3,c), 
G(E) shows a significant decrease throughout the gap for 
strong disorder C5 and C7. 

In summary, a general trend emerges from the ab initio 
two-probe device simulation concerning the robustness of 
the surface helical states of TI material Bi2Se3. For the 
Bi2Se3 2D film with six quintuple thickness, a perturb- 
ing potential of ~ O.leV or less in the scattering region 
of one primitive cell wide, is difficult to destroy the heli- 
cal states. Larger perturbing potential may destroy the 
helical states if applied to regions containing roughly 10 
atoms in each unit cell. A single impurity in each unit 
cell - corresponding to 1.67% density, does not derail he- 
lical states up to SV w leV. The helical states at higher 
energies are much harder to destroy by the disorder po- 
tential. When helical states are destroyed by disorder, 
they are no longer localized on the surfaces but are spa- 
tially mixed throughout the film. Microscopically, when 
propagating helical states pass through a disordered re- 
gion, evanescent states are excited around the perturbed 
atoms and the scattering states are modified. 
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